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mathematical  analysis  of  the  effect  of  the 
lethality  va.  dose  function  on  the  ejected 
number  of  casualties,  and  a  better  ■under¬ 
standing  obtained  of  vhen  certain  approxifflations 
apply.  It  depends  on  the  introduction  of  a 
fandom  ■Variable  for  lethality. 

2)  The  effect  of  the  lethality  ■vs.  dose  function 
on  the  ■variance  of  casimltles  is  obtained.  This 
variance  is  now  e:g;ires8ed  in  terms  of  the  Joint 
distribution  of  men,  the  lethality  distribution, 
and  Joint  distributions  of  the  total  dosages 
Ingested. 

3)  The  possibility  is  considered  of  using  a 
second  Monte  Carlo  procedure  to  compute  the 
A(D)  curve  from  each  sarqple  or  trial  obtained 
by  the  first  or  basic  Monte  Carlo  j^eedure.  A 
series  of  test  eeleidatlons  to  be  eosputed  Is 
given. 
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I.  C— ualty  P>odttotlon 

a)  Tba  Bapected  !hafcer  of  Camialtles  E(C)  « 

In  Appendix  A  ve  outline  the  nathenatleal  end  eonoputlng  etepa 
Involved  in  the  calculation  of  O(x^y)  without  derivation  or  motivation. 
This  la  aort  of  a  matheeiatlcal  flow  chart  for  the  calculation.  The 
varloua  Inputa  and  apsaroxijnatlOna  needed  are  alao  described.  0, 

ve  obtain  E(C)  by  the  Equation 

(l)  E(C)  ■  J'  co(x,y)0(x,y)dx(3y 

x,yeA 

_ If  the  expected  density  of  munitions  ijiu.v)  ia  constant  t] _ thm 

Q(x,y)  la  a  oonstant  Q  and  ve  have 

_-(2) - E(C-)-«-0~-/^-^(D(x>y)dxdy - - 

x,ycA 

In  Appendix  Botes  on  the  Calculation  of  O.I,  we  give  a  series  of 
notes  on  the  details  Involved  in  reducing  the  Indicated  functions,  inte¬ 
grations  and  other  operations  to  numerical  quantities.  E.g.,  the  functional 
Inputs  needed.  CjByoxloaiitloinff 4  possible  assumptions,  alternative  conputing 
technlqiuea,  rooosBNudattoul,  etc..  It  is  ooneoived  that  Appendix  A  gives 
the  basic  natheoatle  structure  of  the  calculation,  while  Appendix  B  gives 
decisions  ahd  details  that  have  to  be  made  In  realizing  this  structure. 
Appendix  A  should  not  require  much  change.  Appendix  B  will  require  ei^paaslon 
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M  non  details  are  added. 

Xa  Appendix  D  «e  deeerlbe  a  netbod  for  laptovliig  or  Blavllfyiag  the 
Mtheawtloal  aoaljals  of  the  affeet  of  lethality  on  S(C) .  fhla  also 
leads  to  aueh  hatter  tmdaarotsBdlBg  of  Shat  oeeurs  ■athSBstleally  stasa 
the  dlstrlbutlOB  of 

(3)  Ihgasted  hy  a  wn  In  If  he  Is  there.  . 

Is  not  lognoraal*  and  henoe  of  vhen  the  assunptlon  of  logaonallty  leads 
to  a  pood  agnroxiastlon  of  0* 


h)  The  Varlaaee  of  Casualties, 

In  Appendix  D,  the  effect  of  lethality  Is  now  included  In  the 
■Btbenatleal  expression  for  0^  .  Xi  previous  work  vas  obtained  la 

terns  of  the  Joint  distribution  of  nan,  and  the  Joint  conditional  distribu¬ 
tion  of  casualties  at  a  point  if  a  aan  Is  there.  Only  mean  and  oowianee 
fhnetlons  are  needed.  The  Inclusion  of  the  effect  of  lethality  Is  obtained 
bgr  expressing  tbs  necessary  mean  and  oovariaaee  functions  of  conditional 
casoaltles  at  teo  points  In  terms  of  the  lethality  ftmetlons  or  distribution, 
and  the  Joint  distribution  of  the  D^. 

If  tb.  U^wtaU  of  «d  U  M«rl«. 

than  the  neoesswry  oaloolntlons  beeome  greatly  simplified.  Baeults  are 
given  for  both  this  anssmptlon  and  tbs  gansral  osss.  Xo  ths  gnMral  esse, 
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a  method  is  given  for  obtaining  a  bivariate  distribution  of  and 

vlth  given  maarglnal  additive  distributions,  particularly  the  gamma  or 
e.d.  The  results  reqiulre  the  calculation,  for  each  1  and  J  of  the  five 

Of  course,  these  five  paraaeters  describing  the  ^Int  dlstrlhutlon 
of  and  Dj  ,  Involve  the  rest  of  the  problem,  l.e.,  the  Joint 

dlstrlbutlona  of  breathing  rate,  the  nmltlon  contaminant  pattern,  and 
the  aunltlon  Impact  points. 

The  slnpllflcatlon  that  occurs  If  t)(u,v)  ,  the  e^qpected  dutslty 

of  sunltlons.  Is  cone  tent  Is  also  given  In  Appendix  D. 

n.  Hants  Carlo  Sleulatlon  of  A(D)  Curves 

tbnte  Carlo  siimil&tlc^  or  eaIculatloh,6f  the  distrlbutlbh  6^^^  eaiualtlM 
or  the  distribution  of  the  A(D)  curve  Is  being  considered.  Ih  such  a 
simulation,  a  field  trial  Is  simulated  on  a  computer  by  selecting  observed 
values,  by  the  use  of  random  nuniben^from  the  appropriate  probability 
distributions .  Vor  each  such  field  trial  ,  the  problem  arises  of  computing 
the  obtained  A(d)  .  This  can  be  done  by  systematic  mathematical  procedures, 
but  -those  tried  to  dstte  require  much  calculation.  The  altstmatlve  possibility 
of  using  a  seoond  Monts  Carlo  procedure  to  cospute  each  field  trial  A(D) 
eozvo  is  under  oonsidsratlon  and  the  vork  is  deserlbed  in  Appendix  C  . 
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APPENDIX  A 

Memo  ^26.03  M.  L.  Norden  December,  I96I 

An  Outlime  of  the  Calculation  of  G 

I  Introduction 

In  the  Fourth  and  Fifth  Quarterly  Progress  Report  (Progress 
Reports  526. I7  and  526. I8)  the  expected  number  of  casualties,  E(c), 
vas  mathematically  derived  In  terms  of  various  known  functions  and 
parameters.  The  motivations  and  derivations  are  given  In  these 
progress  reports.  Me  now  wish  to  work  toward  the  numerical  calcu¬ 
lation  of  E(c)  .  In  this  memo  we  shall  outline  the  mathematical 
and  computing  steps  involved,  (without  derivation  or  motivation)  euid 
describe  the  various  Inputs  and  approximations  Involved. 


_ In_general_we_wl 6h_t  o_comput e _ 

(1)  E(C)  =  J  /(a(x,y)dxdy 

x,y€  A 

where  w(x,y),  the  expected  density  of  men,  is  given  and  G(x,y), 
the  probability  that  a  man  at  (x,y}  would  become  a  casualty,  depends 
In  a  complex  way  on  the  remainder  of  the  problem. 

In  this  memprandiw,  we  shall  assume  that  the  mimltlon  tensll^ 
ij  (u,v)  Is  constcmt.  This  Is,  we  shcdl  assume 

W  tr  (u,r)  *  tf 


5 


3 


This  assun^tlon  makes  a(x,y)  =  a  =»  a  constant.  Hence 

(3)  E(C)  «  a  f  J  a)(x,y)dxdy 

x^ycA 

and  we  now  wish  to  describe  the  calculation  of  G  .  At  a  latter 
time  it  will  be  desirable  to  not  assume  (2),  Then  Q(x,y)  will 
have  to  be  separately  calculated  at  the  different  points  (x,y)j 
but  the  steps  will  be  intrinsically  the  same  as  those  given  in 
this  memorandxim. 

In  this  memorandum,  we  shall  outline  the  mathematical  and 
coo^utlng  steps  and  details  involved  in  the  calculation  of  Q,  so 
“tfiat“^e  choices  involved  in  programming  the  numerical  calculation 
of  G  C6U1  be  examined,  and  the  nece8saz*y  flexibility  of  possible 
ir^uts  to  the  calculation  determined. 


* 

In  Msioo  526.04  we  shall  give  Notes ,  which  will  contain  more 
details,  e.g,  the  functional  inputs  indicated,  approximations, 
assunqitlons ,  alternative  computing  techniques,  etc* 

II  The  Calculation  of  Q 
0  is  defined  as 

(4)  0  -  /  L(D)f(D)dD  -  7  Il-P(D)]  dD 

0  0 


*  Appended  ae  Appendix  B 
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Mfjtiere 


(5) 


L(D) 


1 


e 


-yV2 


dy  =  Log  Dosage  vs  Problt 
lethality 


(6) 


P(D)  =  Probability  that  the  dosage  at  a  point  (x,y)  is 
<  D 

f(D)  =  Probability  density  = 

dD 


The  probability  distribution  P(D)  will  be  mathematically 
specified  in  terms  of  its  mean  and  variance  .  (Notes  A, 

Bj  C,D,E) . 

We  now  assume  the  special  breathing  rate  function  model  of 
the  progress  reports. 

The  quantities  (Ajj  and  are  then  given  by 

m  - - - 

^  t=0 


(8)  M«  +  og 


+  Jig) 


OQ 


/ 

t=0 


t 

/  [1-H(t)|  [Ti®R(t)R(s)+WR(t,s):!d8dt 

8=0 


where  dj^,  a^,  ^  and  u  are  known  constants  and 


u 


if  the  munitions  are  placed  at  random, 
if  the  munitions  are  at  fixed  points. 


(10)  H(t)  =  Prob  [T  <  t)  . 


H(t)  will  be  a  known  probability  distribution,  possibly  a 

thres  parasMter  log  nonsal.  (Notes  P,a,H)  . 
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I 


4 


and  finally 


(16) 


XQo* 


(»-5*l*  +  _ l! _ 

+  V’St  (*“  +  -vN 


Where  X,QQ,u,a.c^,ay,a2,k^,I^,k2  are  given  constanta. 

_ in-thlB_woz4c,_.we__have_le_t _ .t.»a.»x»y«JUL»v _ range  syiabollcally  to _ 

«•  or  e  •».  It  nay  be  neceseary  to  truncate.  NUmerleally  notleable 
effect  a  on  0  should^phyBlcallynot-arlse  from  large  valuea  of 
these  variables. 

fbe  steps  do  not  have  to  perfomed  In  this  order.  An  alternative 
sequence  Is  the  following. 

OegMaiag  the  e taps »  we  ebtain  fer  B(t) 

'4^  t 

(17)  *(t)  -  /  /  /  x(u*v,T)dQ(t-T)dvdu 

ymmm  tM) 


ovdtv  Int^gratloiif  w#  obtain 


^  44i 

(18)  *(t)  -  /  /  /  x(«»v,T)dvdu  dft(t-r) 

u«-** 


-  /  K(T)4ft(t.T)  - 

t«^ 


/  K(t-T)4ft(T) 


•ubstitutlng  Equfttioii  (4),  wtt  obtain 


(19) 


»(t) 


(1-a)  K(t)  ab 


/  K(t-T)#"'*‘^  dT 

imO 


.  (l-a)  Kt)  +  «*»•*’*’ 


/  K(T)«‘*''dT 

‘M) 


(80)  K(t)  -  /  /  x(u#v»t)dudv  -  ■— — 

u—  V^(bP^a,) 
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3.T 


77^ 


dT  ] 


tr  8 

/  If  X(u,V,T,)x(u,V,T2)dq(s 

ymrnm  Tj»0 

*  dQ(t-T^)  dvdu 

?  E(t,  ,T2)dq(»-T2)dft(t-T,  ) 

T^O 


(25)  K(t,8)  -  /  /  X(u*v,t)  X(u#T>»)  dvdu 

Ua^«>  Vm-rnm 


.  gvCtzglV- 

_ 


Aialn«  btt08U8«  of  tho  juap  in  Q(t)  at  t«0>  m  hava  from 


Bquation  (22)  and  (4) 


(a^) 


R(t,8)  -  (l-a)%(t,e)  -t-  (1-a)  ab  a“”  /  K(T^,8)e  ^  dT^ 


Tj-O 


+  (1-a)  ab  0“®*  /  K(t,T«)  a  ^  dt- 

*  ‘ 


♦  (ab)*  /  /  K(Tj,Tg)  •  *  ^  *’4»2,dT 

Ti*0  Tg-O  ^  ^  ^ 


AHMIBB  B  4.1 

Nbbo  326.04  II.L.  Horten  BMOteor^  1961 

MOtos  on  tho  Cnlonlatlon  of  0.  i 

te  Ajpnflte  A,  tf«oo« 


k,2 


(1-3)  Vg 

Nott  C  t  Simplification  If  P(D)  Is  Lognormal 

He  :  Equations  (4)  and  (6)  and  Notes  A  and  B. 

If  F(D)  is  log  normal,  then  Equation  (4)  becomes 

O.-S^  "f  ,-V2  a, 
y/Sr 

where 


M 

and  we  have~ellmlnated'~one~  troublesome- IntegratlonT - 

Hote  D:  Ihe  Relationship  between  and  More  Usual  Problt 

Lethality  Notation 

_  . . .  . . He  Equation  (  5)  . . . 

In  the  usual  Problt  terminology 

Pnsblt  ■  a  *  ¥  IMt  -  5  ♦  ) 

In  this  notation 

l»L  ■ 


-x4- 


(1-1) 


Various  suggsstlons  for  Intsgratlng  Equation  (4)  whan  F(D) 
is  not  log  normal  art  glv«!i  in  the  Fouz*th  Qiarterly  Progress  Esport  [i] 
pages  19  to  23.  In  partieular  Technique  2  seems  appealing. 


Present  data  indicates  that  H(T)  can  be  described  as  follows 
-  - i{a>  -  pj(^  -[*  <1^3-  -  - . . . - . - . - . 

(p.l)  T  «  +  T, 

That  is,  the  random  variable  T  is  the  sxim  of  the 
alarm  time  and  the  time  to  stop  intake  of  the  contaminant.  Then 
if  we  let 

(P-2)  Tp  -  reaction  time  -  time  to  the  stopping  of  breathing 
(MIL  i;,  m  Basking  tte 

it  was  found  that,  particularly  at  hl|^  work  rates,  BaililBg  can  occur 
before  breathii^  stops.  BMioe 
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4.4 


(F.4)  Tg  -  min  [T^  > 

It  was  also  found  that  the  probability  dlstrlbutlm  of 
can  be  approximated  by  a  three  parameter  log  normal  distribution, 
the  three  parameters  being  a,  and  Og  .  That  is,  the  random 
variable  -i-  a  is  distributed  in  a  log  nonsal  distribution  with 

paramsters . 


Ug  “  B  r  +  a)  ] 

(f-5) 

aj  ■  o*  [in(Tp  +  a)  ] 

Yes ,  this  allows  the  random  variable  to  take  on  negatlva 
values.  This  Is  meant  to  correspond  to  the  alarm  being  sounded  while 
exhaling. 

The  probability  distribution  of  Is  also  approximated 
by  a  log  normal  distribution.  Here  the  choices  seem  to  be 


(a)  A  two  parameter  log  normal  dlstx*lbutlon 

(b)  A  three  parameter  log  noxnal  distribution 

(o)  A  three  parameter  log  normal  distribution  with  the 
pavamstsr  "a"  having  the  sasw  v^ue  aa  for  the 
distribution  of  T^,  . 


•HS'I 


t  The  calculation  of  H(T 


In  Not*  P  m  defintd 


1^.5 


(0-1)  T,  -  nln  [Ty  >  i;,] 

nie  probability  of  T.  can  be  eonputed  In  two  ways.  Ihe 
flrat  or  exact  solution  Is  as  follows.  Let  I^(t),  I^(t)  and  ^(t) 
be  the  probability  distributions  of  and  respeotiwely.  Ulan 

(0-2)  1  -  Hg(t)  -  (1  -  Hp(t)3  [1  -  \(t)3 

and  so  the  correct  answer  can  be  obtained.  Note  that  if  I^(t)  and 

are  log  nowAi  .  exactly  los-nornal-. _ 

However,  we  do  feel  that  masldng  rarely  occurs  before  breathing 
stops  (the  exact  probabilities  can  be  obtained  from  the  data),  and 
hence  should  have  on  a  slight  effect  on  I^(t).  So  pexhaps  sane 
approximation  can  be  given  which  is  easier  to  compute.  Hence  the 
second  solution  would  be  to  approximate  the  probability  distribution 
of  presumably  by  a  two  or  tluoe  parameter  log  normal  dlst^^ 

Hast  is  we  would  try  to  approximate 

ii-«i 

or 

(®-4)  I  fin  (T,  +  a))  Mi  o*Un  (T,  +  a)) 
faom  tiM  fiPMiteri  of  T,  and  T,  . 
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Wt  could  th«i  u«o  these  values  to  sf>proxlaate  a  log  nonaal 
dlstrllmtlon.  Such  a  technique  has  been  derived. 

Detezvlaation  of  whether  the  approximation  or  the  exact 
solution  are  preferable  for  calculation  needs  to  be  made. 

Mote  H  t  Be  tBouatlons  (7).  (8),  (9)  and  Motes  P  and  Q 

A  decision  of  what  to  do  about  the  probability  distribution 
of  alarm  time  still  needs  to  be  made.  If  It  Is  assumed  to  be 
a  constant,  this  value  needs  to  be  entered  as  a  parameter. 


Mote  I  t 


t  Controls 


4/r 


4,  o  a  0,  (Tm  °  0  (eonstant  breathing  function) 

5.  ■"  0  (Negleot  maslclng  time) 

Drop  of  Munltlona 


1.  1}  »  TJ 

2.  o  =•  0 


Random  Drop 

Fixed,  known  Intact  points 
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APHBiDn  C 


Nhb  $26.05 


n.l.n.A*>h. 


Jaouwy  19^ 


Maett  Carlo  Bitlaafeloa  of  th>  Aft))  corv.  I 


9>la  BHD  Is  a  preliminary  report  on  some  of  the  thinking  of  inu  and  CSIDL 
personnel  on  the  perohlera  of  determining  the  distribution  of  the  A(D) 
curve  by  Monte  Carlo  teehnl^xes,  vbeare  A(D)  Is  the  area  covered  by  a 
given  dose  or  dosage. 

Conceptually  ve  have  a  field  trial  where  minltlons  lopact^  the  contaminant 
spreads  from  the  nmltions.  It  Is  ingested  (say  by  breathing)  and  a  dose 
Ingested  occurs  at  each  point  of  the  target.  Hence,  A(D)  ,  the  axea 
-oovered  by  a  dose  1)5  can  be-obtalned- for  the  trial.  If  “another'  trlgl''i¥~ 
zun  under  "Identical  conditions"  the  munitions  nay  Impact  at  different 
points,  etc.,  and  another  A(D)  curve  will  be  obtained.  If  very  many 
trials  are  run,  at  least  conceptually,  ve  will  obtain  many  trial  A(d) 
curveiV  lahd"^  can  be  describe  by  a  prob^lll'^  distribution  of 
A(D)  . 

Tor  given  mathematical  pDrobablUty  diatrlbutlons  or  exact  deacrlptlona  of 
munition  Impaota,  contaminant  dlsaemlnatlon  patterns,  breathing  rates 
(and  positions  of  msn  and  iethailty^^^f^  if  ve  vent  casualties):,  ve 

can  proeesd  In  two  vays. 

Itoft  a  MthaBstlnal  daseriptlon  of  tbs  prohsiblllty  dtstrlhutlM  of  A(]>) 
can  he  Obtained,  and  voA  In  proceeding  along  these  Unas. 


•  If  Ingsetlon  Is  negleeted  ve  have  dosage.  If  the  dlstrlhotlon  of  aen 
and  a  letlinllty  distribution  are  added  as  li^uts,  ve  vlU  have 
oasualtlse. 
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Ihe  seeand  motbod  !■  to  nodal  tba  trials^  usliig  Monte  Carlo  prooedurea, 
on  a  computer,  hence  obtaining  the  e^il^railent  of  a  nuoiber  of  trlale  and 
hence  of  a  nuMber  of  trial  valuea  of  A(D)  ,  vhleh  can  then  be  uaed 
to  eatloate  the  probablUtgr  dlatrlbutlon  of  A(l))  .  Thla  la  done  bgr 

aeleeting  random  nunibera  from  the  appropriate  probability  dlatrlbutlcm* 
Tor  example,  if  aeven  nunltlona  are  dropped  on  one  aiming  point  vlth 
a  known  probability  dlatrlbution  of  laipaeti  then  ve  can  seleet  aeven 
random  lavaet  polnta  from  this  dlatrubutlon  and  uae  theae  aeven  polnta 
in  a  computer  trial.  Similarly  all  other  valuea  for  the  trial  are 
aelaeted  randomly,  and  ve  obtain  a  mathamatiecCJy  defined  and  coaputable 
A(D)  curve  for  the  trial.  Thla  proceaa  la  repeated  until  a 
aufflclent  nuid>er  of  triad.  A(D)  eurvea  are  obtained.  Let  ua  call 
thla  the  flrat  (or  baalc)  Monte  Carlo  paroeedure. 

At  this  point,  the  queatlon  arlaea  of  how  to  eonpute  the  A(d)  curve 
maaerandum. 

Thla  calculation  haa  been  tiled  by  a  ayatematlc  procedure.  The 
poiaibllity  of  \iiing  anothof  Monte  Carlo  procedure  to  obtain  the  A(d) 
for  a  given  trial,  la  nov  being  ezaadned,  a  aort  of  aecond  Mnxte  Carlo 
procedure  vlthln  a  flrat  Monte  Carlo  procedure. 

The  ayatematlc  procedure  used  breaks  tha  target  rectangle  up  into 
auRll  areaa  of  aiee  dacdy  and  caaq^utes  A(d)  by  counting  the  nunibar  of 
au^  araaa  vlth  a  central  dose  >  D  •  Here  the  error  is  a 
"nathaeatleal  error*  and  depends  on  doc  and  dy  .  jy  "error*  ve  mean 
tba  dlffiKeniea  batwein  tha  eorreot  A(d)  curve  and  the  A(D)  einvl 
calculated.  fOr  given  numarleal  valuea  of  D,  dx,  dy  and  the  atartlag 
points  of  the  Interval,  the  error  la  a  eorreapoiiding  nuMber,  asy  e(D)  . 
Ve  eall  thla  naiamr,  a(D)  the  "mathamatioel  anor*.  Howvar,  for  tbls 
proeedure,  the  nuMbar  of  caleulatlona  also  dapanda  Inveraely'on  dx 
Md  4r  •  va  find  that  If  ds  and  dy  are  large,  the  maiiar  of 
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caleulatlons  xieeded  can  be  aaell,  but  the  "nathOHatleal  error"  vlU  be 
large.  SlBdlarly«  If  ^  and  Ay-  are  nade  aaall,  the  aetha—tlcal 
error  can  be  made  nallf  but  at  the  expense  of  a  large  number  of 
■Bchlne  ealeulatlona. 


One  alternative  la  to  try  to  develop  other  vay  of  .eystematlc  nutaarleal 
ealoulatlon  of  A(D)  .  HSm  ve  ean  eltheTj  (a)  .try  to  lispaEVve  the 
accuracy  of  the  above  procedure  or,  (b)  attaoqpt  to  develcfp  eifentlally 
different  ayateaMtlc  procedure.  As  an  ezaaple  of  (a)  ,  the  above 
procedure  coa^utea  D(x^«  y^}  at  a  set  of  points  •  Ida,  y^  <■ 


If  ee  look  at  the  suberees  of  ehleh  these  are  the  centers,  it  la 
clear  that  aama  Iqprovegwnt  around  D  [irhleh  oecors  at  (O.O)]  adght 

BBiZ 


be  obtained  by  using  the  points  •  u  +  y^  ■  v  Ih 

particular  the  values  u  ■  v  ■  o  or  u  ■  ^  "  ^7/2 

Indicated.  As  an  erf’p't*  (b)  ,„we_iBl^t-jconalder_tr7lng_tOLjOGBipute 
or  approxlaate  the  contour  in  the  (x,y)  plane  eorrespondlng  to  a 
given  D,  and  then  nuMrleaUy  Integrate  to  obtain  the  area  vlthln  It. 
This  has  been  tried  using  slngple  Interpolation  to  obtain  -the  contour, 
-and  gavs-a-largar  ermr.  Posaibly-the-use-of-higher  order  Inter¬ 
polation,  or  the  of  appropriate  curves  to  the  oontoura  (say 

hlvarlate  polynoaleela)  viU  give  better  results. 


A  seoond  eltematlve  is  to  use  a  second  Mcmte  Carlo  procedure  to 
ooavute  A(d)  .  ■ere  the  nui^er  of  eelculations  n  and  the  siae  of 
doc  and  dy  are  unrela-tad.  n  Is  usually  called  the  seavle  siae, 
■id  for  any  dx  and  dy  ee  can  seleet  any  sas^le  else  n  ve  vtah. 
■seever^  ve  nor  hqnre  a  ’^■thaaatleal  errcr*  dspenrtlng  on  ds  sod 
and  a  stattstloal  tnrar  depsndlng  on  r  •  Ihst  la,  tbs  ssooad  Manta 
Ovlo  proeadure  resolta  In  an  a(D)  cures  eceputed  trm  the  n  values 
of  D  Qibteined.  Boeevav  unlike  the  systwntin  pneoedure,  if  the 
•aonnd  Mnte  Carlo  proaens  er  eaToulatlon  is  repeated  on  a  Bonpitsr» 
a  dlffOrat  set  of  rsniSB  nuMbsrs  vUl  he  selaoted,  henee  a  dlfitaront 
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■at  of  n  yaiLues  of  D,  wd  hanea  a  dlffarant  A(D}  eurva  vUl  te 
oibtainad.  That  IM,  tha  A(D)  obtalaad  ia  vow  a  randoai  nerlaibla.  fbr 
any  D,  A(l>)  ,  vhlla  banriag  an  avarasa  or  aapaatafl  aaljoa  I  [  A(D)  ] 
and  a  irohahlllty  dlatrlhutlon  aibout  S[  A(D)  ]  vhleh  ean  ha  daaerlhad 
In  taxna  of  o^[  A(D)  ],  tha  ataodaEtd  derlatlan  of  A(D)  .  B[  A(D)  ] 
1»,  for  axQT  obDlee  of  4ac  and  dgr  ,  tha  ralna  that  aould  ha  ealaulated 
hf  tha  aibova  ayataantlo  gacooadura  haaad  on  tha  aana  dn  and  • 

Hanea  tha  dlfferanea  hataaan  B(  A(D)  ]  and  tha  coarraot  ralua  of  A(D) 
la  the  anthegnttleal  error  daaerlhad  aibova  (eallad  tha  '^aa"  In 
■tatlatleal  tandnoloar).  Mote  that  It  doaa  not  dapaad  on  a.  Tha 
"atatlatleal"  error  of  tha  A(D)  ealeulatad  la  daaerlhad  in  taraa  of 

o^t  A(D)  ]  and  dependa  On  n,  not  on  dr  and  dgr  • 

In  faet^  alnea  the  randOB  numbera  uaad  to  aalaet  the  polnta  (z^y)  ara 
al£d>t  digit  decimal  a  >  It  turns  out  that  It  aaor  ha  noat  eonvanlant  to  taaa 
a  ^  and  dgr  of  O.OOOOCXX)!,  vhleh  la  vary  anall  and  rasulta  In  vhat 

error  of  tha  saoond  Monte  Carlo  proeadure  la  Juat  tha  atatlatloal  one. 

It  should  ha  pointed  out  that  even  for  tha  saaaa  Impact  polnta,  ate*, 

~tha~  saeond  Monte  Carlo  xarooedura  (aaaiipllnB  of  tfhaarvatlon  polnta)  vtll 
Introduea  varlablU^  In  the  A(D)  curve*  ■owvrai,  for  angr  eat  of 
nhnaan  d's  thla  ean  ha  reduced  by  aanpUng  aora  polnta.  The  aaln 
lonhlan  la  to  datanalna  Juat  hov  nany  polnta  to  aalaet. 

It  la  danaad:  prefaraibla  to  eonoantrata  on  the  aeoond  part  of  the  prohlaai 
(aatinating  the  A(D)  curve  for  given  drop  data)  before  going  to  tha 
flret  part  (aatinating  average  A(D)  curve  and  Its  variaiblllty  for  a 
l^van  taatiaal  ^tuattion)  *  Xn  tha  Taiigmp  cf  Iteflihaitte  pp 
thla  aaaood  prohlan  involvaa  treating  a  alngla  realiaatlOB  of  tha 
doaaga  field,  rather  than  tha  ansaaibla  of  flalda*  The  varlaibllity 
Obtained  In  thla  pliaaa  will  have  to  ha  added  to  the  variability  dna 
to  Inipaet  polnta,  ato.  (first  phaaa  of  the  prohlan). 
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oplortt  tbs  posBlblUtlas  of  using  such  a  ssoond  tgrpe  of  Monts 
Carlo  proesdurs,  tbs  foUovlBg  nsehlns  oaloulctlOBS  vsrs  agrssd 
tqpon. 


1.  lbs  slJMdo  mition  asttsm 


(1) 


fiks  tfsa  eovarsd  by  D/q  -  0  Is  to  bs  obtalasd, 
D  2  cr  X  ^ 


And 

Cy  -  0.211  Mtsrs 

C.  «  0.135  ostsrs 
s 

u~^2~»tiw[/i8beea  -  -  — 

a  •  ^  ■  0.25 

S  is  la  xuilts  of  eem  ssoonds/Cnstar)^ 
Q  la  la  ualts  Of  gras* 


Itaoe 

0  Is  la  ualts  of  sseonds/Castsr)^ 

A  total  arsa  0»lag  fiwB  x  ^  O  to  250  BBtsTS  la  tbs  x  ooordinata 
and  froB  -16  to  ’flB  ■stars  la  tbs  j  ooordlaata  la  oonsldarad 
[  9,000  sq.  ■Btapsl  . 


■ota:  fba  axaot  A(d)  dorva  ftr  tbls  pattara  e«i  ba  obtalaad 
■atbsaatlealiy  and  will  ba  usad  fbr  oo^^arlson  vltb  tbs 
BUMTlaal  osleolatlons  to  ba  dssarlbad. 
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2.1.  aad  lover  dosetK  ^toffv  viU  le  used.  Tbm  deitdlle 


2.2  fbft  tidiMt  of  B-  are  then  tabled  Istb  100  olaM  latarvali 
giving  A(9)  at  100  valnea  of  B.  The  Mlaetion  of  theae  elaM 
Intervala  la  being  enaadnad  In  nore  detail.  One  point  la  that  the 
loveat  interval  aee«n  to  define  a  aaeond  kind  of  aero**. 


*  Otfer  and  lover  doaaga  eutoffa  iuat  be  uaed.  Equation  (l)  InpUea 
that  doaage  la  poaltlve  for  all  z  and  y  and  can  eaally  give  valaea 
greater  than  the  ■aehlne  can  handle.  Doaagaa  above  a  lethal  value 
een  be  reduced  to  the  lethal  value  without  any  Iona.  It  la,  of 
oourae,  alao  neceaaary  to  aiit  down  the  alee  of  the  doaage  field 
ao  that  it  can  fit  In  the  counter.  Thla  baa  the  effect  of 
truncating  at  the  low  doaage  yaluea. 

**  The  first  kind  of  sero  occurs  vImk  the  doaage  field  la  "bozed**  In 
by  a  finite  "target”.  Thla  aecond  kind  oceura  when  values  are  left 
In  because  of  the  shape  of  the  contours.  See  Figure  1,  vheare  we 
—avwoaa  that-  0  M  .00$-is-the--jrtni»ai-iaDjBtour.  — T^-1^ 

the  reetflOBi^  repnieent  the  first  kind  of  eero,  and  x  rapreeenta 
the  second  kind,  which  oonalats  of  the  points  between  the  bounding 
reetengla  and  the  winlviM  oootour. 
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A  second  point  Is  that  logarithmic  intsrva;ls  in  6  would  he  preferahle. 

For  example,  lOO  logarithmic  intervals  from  3  to  ^00  would,  make  each 

division  about  la047  times  previous  division.  If  we  make  the  ratio 

1.0^  (that  is  5^  increments)  then  100  logarithmic  intervals  gives  a 

# 

ratio  of  about  132  from,  the  fiarst  to  the  100th  Interval  . 

3o  CauLculationa  of  A(fl)  by  use  of  a  systematic  grid  have  been  per“ 
formed.  It  was  decided  to 

3.1  Repeat  the  systematic  calculations  for  the  single  munition 
for  the  two  subareas 

a.)  =  50  ,  ^  =  2  and 

b)  /:Uc  =  1  ,  =  .5  o 

3.2  Perform  Monte  Carlo  calculations  for  A(e)  for  a  single 
munition  as  foU.ows* 

3.2.1  Eight  digit  random  numbers  were  to  be  used  making 

_ _ _ _ -/V“^-tOOOOOooi---; - -  ^  ^ 

3.2.2  &-\00  hundred  points  were  to  be  conqputed  in  the  follovrlng 
manner.  16  runs  of  400  points  each  were  to  be  made,  and 

_ _ jtoe  resulting  A(_0)  tabulated  in  the  lQ0„_clasB  intervals 

for  each  run.  The  successive  4  runs  were  to  be  tabulated, 
giving  4  runs  of  l600  points  each.  The  total,  1  run  of 
6400  points,  is  also  to  be  tabulated. 

3.2.3  Each  of  these  estimated  A(e)  (16  from  400  point,  4  from 

. -  1600  polnts>  1  from  6400  points)  la  to  be  compared, 

presumably  graphically,  with  the  exact  A(0)  which  is 
known  In  the  single  munition  situation. 


*  Since  d  a  D/q,  corresponding  changes  in  log  9  are  the  same  as 
corresponding  changes  in  log  D.  However,  o  ehoBoe  in  D  is  Q  times 
the  corresponding  change  in  9. 
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5.2.4  TOie  curve*  A(0)  +  c[A(e)]  and  A(0)  +  2c[A(9)]  are 
also  to  be  ccniputed  and  presented  fCr  conqparleon* 

3.3  Perform  tbnte  Carlo  calculations  similar  to  tbpM  made  In 
Section  3.2  for  the  two  flmad  aunltlons  situation. 
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MeBD  526  o06  MoLcN.  February  I962 

The  Concept  of  a  Random  Variable  for  Lethality^  Its  Bffect  on  E(c)  , 

p 

and  the  Inclusion  of  the  Effect  of  Lethality  on  a  (C) 

1.  Introduction 

In  Progi'ess  Reporta  526ul7,  5260I8,  526,19  [l]  [2l  [3],  a 
mathematical  model  of  the  structure  of  casualty  production  was  given, 
the  expected  number  of  casualties  E(C)  derived,  and  \rork  toward 

c 

the  variance  of  casualties  described.  In  this  memo  we  shall  describe 
further  improvements  in  the  mathematical  analysis  of  the  effect  of 
lethality  on  E(C)  and  continue  the  work  of  [3]  to  Include  the  effect 
of  lethality  on  . 


In  Eqmatlon  (25)  of  Progress  Report  526. 17»  the  probit  versus  log 
dose  lethality  function  was  defined  as 


InD 

(1)  L(D)  -7=  f  exp  {-(x-<*v)^/(2(Tj^^))  dx 


1 


L 


-y% 


dy 


where 
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(2)  L(d)  =  ptrobablllty  that  a  man  vho  has  Ingested  a  total  dose  D 


will  he  a  casualty. 

and 

(3)  =  logarithm  of  the  50^6  lethal  dose. 

(4)  CL  =  the  reciprocal  of  the  slope  of  the  ]'j(d)  versus  D  curve 

li 

when  plotted  on  log  normal  or  log  prohlt  prohabillty  paper. 

Then,  in  Equation  (3)  of  Progress  Report  526. 19,  we  defined  the 
rand<^7yariable^ - - - — 

(0  if  a  man,  assuming  he  in  A.,  is  not  a  casualty 

(5)  ; 

^1  if  a  man,  assuming  he  is  in  A^,  is  a  casualty 
Then  we  defined 


(6)  0, 


and 


(7){  F.(e)  ■  pirohahlllty  that  the  random  variable  D.  is  <  the 
1  number  6,  *” 

U(«) .  ^ 

^  d6 

♦  The  defined  in  equation  (13^  and  used  on  page  11  of  Ingress 

526.17  is  not  the  same  but  rather  corresponds  to  the  conditioned  random 
variable  given  that  the  dose  ingested  in  A^  is  D^,  The  used 

on  page  12  of  526*17  is  Uie  unoonditlonal  random  variable,  as  we  now 
define  . 
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Then  on  pages  11  and  12  of  526.17,  we  showed  that 

CO  CO 

;8)  Gj  =  J  L(D)f^(D)dD  =  J^L(D)dFj(D) 

0  0 

00  00 
^  J  [1-F^(D)1  dL(D)  =  J  dD 

0  0 

end  pages  15  to  24  are  devoted  to  a  discussion  of  methods  of 
conpitlng  G^. 

II,  The  Introduction  of  a  Random  Variahle,  h^,  for  Lethality 

The  Introduction  of  a  random  variable  for  the  lethality  consider^ 
ably  simplifies  the  proof  of  some  of  the  results,  and  leads  to  the  next 
step  after  P.R.  526.19,  the  inclusion  of  lethality  into  the  calculation 
of  0  The  idea  is  as  follows.  Let 

(9)  h^  s  a  random  variable  called  lethality. 

That  is,  h^  is  a  random  variable  with  a  cumulative  probability  diatrlbu^ 
tlon  function  L(D)  or  L(@).  That  is,  ware 

(10)  Prob  [h^  <  e]  -  L(e) 

*  Again  a  alight  ootatlonal  inconslatency  arlsea.  The  symbol  h^  was  used  in 
PJi.  526.17  fbr  the  constant  a^  to  be  defined  in  Expiation  (18)  of  this  report. 
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e.k 

and  the  h^  sure  statistically  Independant  of  each  other  and  of  all  the 

other  random  variables  in  the  casualty  model. 

Then,  if  we  define  as 

(11)  kj  = 

the  model  is  Just  as  before.  We  have  just  intraduced  a  symbolism  for  the 
lethality  random  variable  rather  than  just  its  probability  distribution 
L(0)  . 

III.  The  Use  of  Tills  Nei;  Random  Variable  for  the  Calculation  of 
The  conditional  expectation  of  given  is 

(12)  E(kjlD^)  =  Frob  (k^  «  iJd^)  =  Prob  (h^  <  |d^}<^L(Dj) 

and  the  unconditional  expectation  is 

CO 

(13)  =  E(k^)  =  Ef.L(Dj)3  =  J  L(D)dP^(D) 
aa  given  by  Equation  (8) . 

However,  thla  Ibnulatlon  slao  glvea  ua  other  resulta.  We  can  write 


0  if  h.  >  D. 

1  1 

1  if  h^  < 
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(ll^)  -  E(lc^)  =  Troh  (lij_  <  D^)  =  JJ  P^(d)  dDdh 

h‘•^D 


Now  the  diatrlbution  L(S)  of  h^  was  assumed  to  be  probit  vso  log 

dose  or  wliat  is  commonly  called,  log  normal »  That  is,  the  distribution 
of  Inh^  will  be  a  normal  distribution  c  In  pea-ticulax,  from  Equation  (l) 

(15)  E(lnh^)  =  lij. 

a^(]jihj^)  =  cf^ 

This  also  ij35>lles  that  h^  is  non  negative,  that  is,  it  lies  in  the 
range  0  -bo  +  «o  D^,  by  iliE  veiy  ‘definltloh,  l3 ”1^ 

Hence  we  can  also  write 

(16)  ==  E(k^)  =  Prob  (h^  <  (mh^  <  InD^^) 

=*  Pirob  (0  <  InDj^  Inh^^]  =»  Prob  [0  <  z^] 


where  we  have  defined 
(IT)  -  IbD^  -  Inh^ 


Bov  let 
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distributed  and  even  to  understand  when  Equation  (20)  Is  a  good 
ajqpraxlnat  Ion . 

V.  Calculation  of  when  is  not  Log  Normal,  and  the  Effect  of 

Aroroxiinatlng  by  a  Log  Normal  Distribution 

If  the  probability  distribution  of  Is  not  lognormal,  we  can 

obtain  the  distribution  of  as  the  distribution  of  the  difference 

between  InDj^  and  Inh^^,  and  then  Integrate  to  obtain  by  Equation  (l6)o 

The  result  can  be  placed  In  the  form  given  by  Equation  (45)  of  P.R,  526,17^^^, 
Another  approach  is  described  in  Technique  5  on  page  23  of  p,R,  526. 17. 
That  is  to  errand  the  distribution  of  InL^'  in  a  Gram  Charller  Series  <, 

Here  there  are  two  points.  The  first,  as  pointed  out  in  526.17,  is  that 
we  feel  that  the  distribution  of  InD^  is  "approximately''  normal.  Surely 

it  Is  -HBich-"eloeer*'  to-a -nonnal  distribution  than-  is-  the  distribution  of 

Df. 

However,  ve  now  see  a  second  point.  Since  z^  Is  the  difference 

between  two  random  va^ables,  one  ot  vUch  is  normally  distributed,  the 
distribution  of  z^  will  be  such  "closer"  to  a  normal  distribution  than  is 

the  distribution  of  InD^.  We  can  see  this  by  Io<^l£lxi£  at  the  r'th  cumulanti 

of  the  mdiOB  variables. 

For  a  mmally  distributed  randon  variable,  all  the  cumulants  above 
the  second  are  aaro.  That  la 


'35- 
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(21) 


Kg  -  a‘= 


K  =  0  If  r  >  3 
r  — 


This  defines  a  noraial  distribution. 


Standardized  (or  unitless)  cunulants  o:  are  defined  as 

r 


The  closeness  of  to  zero  (for  r  >-3)  Is  the  appropriate  criteria  for 
closeness  to  nonaallty. 

In  this  notation  we  have 

(23)  Kj^(lnhj)  -  Hj. 

--  l^(lnh^)  .  4  - - 

K^{lnh^)  -  0  If  r  >  3 

(24)  Kj^(liiD^)  - 
KgdaD^)  .  oj 

l^(lnD^)  aa  givan  'bgr  the  dlstrlhutlon 


36- 


6.9 


Hence, 


(25)  aj.(lnDj) 


Kr(lnD^) 


if  r  >  2 


Now,  by  the  additivity  of  cuoulants  of  the  sum  (or  difference)  bete^een 
independent  random  variables,  we  have 


(26)  Kj.(lnz^)  -  Kj.(lnD^  -  Inh^)  =  K^(lnDj)  f 

Hence 

(27)  Kj^dnz^)  =  hi  - 


KgClnZj)  *  +  ej 


Kj.(lnz^)  »  Ky(lnDj^)  «  oJa^(lnD^)  if  r  >  3 
And  so,  if  r  >  3 

"i  T‘'^ 

jr^-J  v".i 


a^dnDj) 

“  (1  •> 

Clearly  Q!  (lnz. )  is  less  than  a  (inD,  ). 
r  i  y  1 


K  (inz. ) 

(28)  a  (to.  )  .  - !_ 

'  ‘  (l^(to.j)f''2 
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The  BmIg  Bgaatlon 

Id  the  last  grogresa  report.  Progress  Beport  ^26.19»  the  vwrlanoe  of 

easualties,  ^  ves  ohtalsed  In  terms  of  the  Joint  distribution  of  nsn  end 
c 

the  Joint  (oondltional)  distribution  of  easualties  at  a  point  (x,y)  if  a 
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man  Is  there.  EssentlcLLly,  only  nean  end  povarlance  functions  of  the  two 
Joint  distributions  are  needed. 

If  the  target  is  broken  into  small  subareas  the  Joint  distribution 
of  men  vas  described  by  the  quantities 


(32)  /  B(k^)  »  ki  -  0^ 


a(kj:  kj)  .  S(kjkj)  -  OjOj 

Ckj-ij  kj"i)  .  *(kjkj) 


A  continuous  representation  was  also  obtained,  and  for  this  the  funetloim 


w(xj^y),  O(Xj^y),  Kjs(*i»yiJ*2>y2)/  ®*«*  Introduce^  where 


«(xj,yj)  - 

0(xj,y^)  -  Oj 

V  , 

-  (rfkjjkj) 


■39- 


For  loth  the  discrete  and  continuous  rejresentatlons  a  was  obtained 

In  tezms  of  the  above  quantities.  A  general  solution  was  obtained  and  a 

nuniber  of  special  cases  given.  At  various  tines  In  the  argument,  the 

derivations  In  tenas  of  the  bivariate  probability  functions  P  and  P, 

n  K 

seem  surer,  and  these  were  used  to  obtain  E(C^).  From  which,  of  course, 
we  could  Obtain  by  the  equation 


i3^)  0^  -  E(C^)  -  CE(C)]^ 


Table  1  lists  these  results. 


table  1 


tlon  Numbers  of  Results  In  Progress  Report  326.1 


aivlng  or  E(C)^  In  terms  of  the  Probability  Distributions  of 


^Hsin  ahd'lJohditional  Cesusiltles  Ohder  Various  Assnniiptlcms 


In  this  vork  the  random  variables  k^  depend  on  the  lethality,  the 

breathing  rate,  the  munition  contaminant  pattern  and  the  nunition  lopaet 
points,  and  hence  susmarize  everything  but  the  placement  of  men. 

The  Joint  distribution  of  the  k^^  enters  into 

mean  0^  or  and  the  covariances  o(kj^:kj)  or  Kj^(xj^,y^jX2,y2) 

equivalent  functions.  The  calculation  of  G  has  been  described  in  the 
previous  progress  reports.  Let  us  nov  obtain  the  covariances  in  terms  of 
the  distribution  of  lethality  h^  and  the  dose  Ingested  D^. 


through  the 


NOW 


(35)  *k^*i*yi^*J»^J^  '  of(ki:kj)  »  -  0^0^ 


■  -  OiOj 


and  ve  see  that  it  will  be  sufficient  to  obtain 


(36)  E(kjkj)  -  Erob  (kj-ljfcj-l)  - 


Ba«r>  ftroB  Bqjuatlon  (U),  ve  have 

(37)  ®(Vj^  "  (V^;kj-1)  -  Proh(D^  >  h^;Dj  >  hj) 


3y  the  noonegatlvlty  of  the  and  ,  and  the  definition  of  , 
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Equation  (l6),  this  becomss 

(38)  E(kjkj)  »  Prob  [D^  -  >  Oj  Dj  -  hj  >  O] 

■  Rrob  >  0;  >  0] 

*  Prob  [z(xj,y^)  >  0;  z(xj,yj)  >0  ] 

These  are  the  basic  equations  ftrom  vhlch  E(k^kj)  Is  to  be  obtained. 


VH.  E(k^k.)  If  the  Bivariate  Distributions  of  D.  and  D.  are 


Now  If  the  Individual  probability  distributions  of  the  are  log 
normal,  then  the  Individual  distributions  of  the  z^^  are  normal.  If  we 


make  the  slight  additional  asBUB^lon  that  InD^  and  InDy  have^a^ 

bivariate  normal  distribution,  then  z^  and  z^  will  edso  be  so  distributed 

(Remember  that  the  h^  were  assumed  to  be  statistically  Independent  of 

each  other  and  the  other  random  variables,  but  the  were  correlated 

with  each  other.  Substituting  Equation  (38)  Into  the  standard  bivariate 

nonil  sorababillly  denUty  flih^ldh,  '%nt  bbt^ 
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(42) 


6  ,i6 


2  2,2 

\  =  “l  ■'  °L 


Now  to  ottaln 


(43) 


r- 


V.1 


Note  that 


(44)  -  *  Cf(lhD.  -  Inh. ;  LnD,  -  liih.) 

.  1  1  j  j 


If—we-defiae- 


(45)  =  eT(lnD^,lnDj) 


Then 


if  1  /  J 


if  1  -  J 


Hence 


if  i  J 


if  1  -  J 


i 
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(49)  -  (inD^jluDj)  s  Oj 

I  - 

Then 


«  i  ^  J 


If  i  =  J 

This  result  Is  very  useful  and  convenient  because  the  function  M(h;k;r) 
has  been  extensively  tabled  and  it  is  believed  that  adequate  ccaqputer 

fMgRMa  miA  afpwi»attcaa  an  aHiliitti  for  ita  ealaulattaii. 

One  of  the  first  tables  of  N(h,k,  )  Is  given  by  Karl  l^eursen  on  pages 
78  •  157  of  [4].  The  aost  extensive  table  Is  given  in  [8]  where  M(h,k,r) 


-45- 


is  called  L(h,k;r).  Extensive  bibliographies  are  given  in  [6]  and 
[8].  A  discussion  of  computing  techniques  is  given  in  C4][^][6][7] 
and  [8]. 


Vnio  A  Simplification  if  n(u,v)  »  the  Expected  Density  of  Minitlons, 
is  Constant. 

O 

If  Tj^  or  Tj(u,v)  is  constant  then  the  calculation  of  or 
E(kj^kj)  is  greatly  slngplifled  in  two  vays;  in  the  calculation  of 
and  in  the  calculation  of  M.  We  can  see  this  as  follovs. 


The  distribution  of  or  D(x,y)  will  now  be  the  sane  at  every 


point  on  the  target.  The  correlation  coefficients 


or  the  co- 


variances  will  now  only  depend  on  the  difference  between,  the  two 

points  (x^,yj),  (Xj,yj),  that  is  on  -  Xj  and  y^^  -  y^  .  Hence  the 

c(kj:kj)  will  have  the  same  property.  That  is,  the  functions  ^ 


and  "  '  . . —  - - - . . .  - 

will  be  functions  of  the  two  quantities  the  four 
qjuoitltles  slavllflss  the  oaleulstlon  of  or 
B(k^kj)  given  in  Table  1. 


Xn  additloii,  if  ve  aasitme  that  the  diatrihutions  oj;’  the  sxe  Miml, 

the  M  calculation  that  now  enters  into  the  calculat-ion  of  each  E(k.k  )  ie 
also  Bimplifled.  That  is,  If  ti(u,v)  Is  constant,  we  have  for  InD^ 


(51) 


for  all  i  ctnd  J 

for  all  i  and  . j 


and  similarly  for 


for  all  i 

X-52-)  — - _  _ _ r 

~  -  ef  +  0?  for  all  1 

Z ,  Z  IJ  i. 


'riien7  if  the  distribution  of  and  z^  aj.'e  bivariate  norital,  we  have 


(53)  =  M(h,b,r) 


where 


(■5V) 


2  2 


and 


(55) 


z 


1  + 


2 


Of  course,  (j.  .  =  rt(lnD.  :lnD.)  and  /O  .  will  be  functions  of  the  'mo 
quantities  and  'the  four  x^,x^.,y^,yj r 

;Now  the  calculation  of  M(h,bjr)  or  M(h,hjr)  is  even  siHX'le^r  tiien 
the  calculation  of  M(h,hj^)  and  considerable  simplification  occurs.. 

For  example,  from  Equation  5  of  f^53>  obtain 


(56)  M(h,h,r)  =  2M(h,0;  « 'fh^r/s) 


Froni  the  Equations  on  pace  7*^  of  [6],  we  obtain,  for  h  >  0 


M(h.h.p)  »  2T(h.'»)  2T(h,  ) 

'f  x4-r 


Where 


(58) 


h^(n-x^) 


T(h,a)  -  ^  /  2. 


dx 


1+x 


For  negatiz-e  h,  we  use  the  result 


48. 
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2 

(59)  ll(-h,-b,r)  -  M(h,h,r)  +  ^  /  e“*  dt 

-h 

These  are  typical  of  the  kinds  of  sisqplifications  available.  However 
the  M  function  can  he  e^onressed  in  other  ways,  and  other  fOrsulations  of 
M(h,h,r)  Bsy  be  preferable  for  calculations. 

Because  of  this  great  sia^lification  available,  ve  strongly  feel  that 
2 

calculations  of  a  should,  at  least  at  first,  be  based  on  the  above 

C>' 

.  bl«rlat.  Jagnorml  dlrtrlbutlon  of  Dj  D^. 

—  l^^he  next-sectira  weUBair ^w~«f8alts“f or “genaal^bivarlat e- 
distributions  of  0^  and  Dy  This  can  then  be  used  to  check  the  degree  of 

approzlBiation  to  o^  obtained  by  the  use  of  the  lognoroal  assunption. 

K.  B(k^kj)  in  the  General  Case,  in  Terns  of  the  Bivariate  Distribution 
of  and  o 

If  the  distributions  of  the  are  not  multivariate  Log  Honaal,  then 

S(k^k^)  Is  iuch  more  difficult  to  obtain.  We  shall  give  a  genere^  result, 

but,  of  course,  viU  be  thinking  of  nunarleal  calculations  baaed  on  the 
gsHa,  e.d.  or  sljsllar  distributions.  Bote  also  that  B(k^kj)  only  depends 

on  the  bivariate  distribution  of  and  Dy  Let  us  introduce  the  following 

notation  for  this  bivariate  distribution. 
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(®>)  Dj  <  SgJ 


<i^4  .(0,;0o) 

(a)  s- 


(62)  ^  ^  ®2^ 


Then  ve  can  obtain  B(kj^kj)  In  two  vaani.  The  first  result  Is 
obtained  by  looking  at  the  oondltlonal  expectation  of  holding 

0^  and  Dj  fixed.  Ey  equation  (37)  and  (lO),  ve  obtain 


(63)  S(ki;kj  |Di»Dj)  • 


Hinde  tlM  imdmid^  eigpeetatldh  la 


W  *0^.1. j)  .  f  f 

.  V»J  . 


»1»2 


6. 


The  •eooBd  remit  is  ohtslned  hgr  first  holdliig  and  flmd.  Than 
the  oondltloiisl  sapsotaitlon  of  k^k.  Is,  hx  Ivuctloii  (37)  and  (6e) 


the  last  step  hy  BqjLiatlon  (^). 


Then  tbe  unoonditloiial  eapeetatlon  la 


_(^)__l(kjk.l 


1 1 


»1»2 


These  two  reaults  depend  on  eaqpreaalng  the  hlyarlate  dlatrlhutlon  of 
and  Dj  In  teaw  of  tbe  funetlona  ^1  ^ 

K.  .(d.  vO  *  These  mat  now  he  obtained  if  these  results  are  to  be 

useful. 
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X.  A  Mrthod  fiar  Obtaining  BlYTlate  DjitrtbiitlOM  of  ani  ghat 
Exn  Olven  Additive  IfcrgliMd  Pl«trlbutlOiiM. 

Bov  if  the  distribution  of  is  asstned  to  he  g— or  e.d.,  there 

is  DO  ual^ie  blvarlste  dlstrlhution  vlth  a  gama  (or  e.d.)  dlstrihatlon 
fOr  the  marginal  dlstrlhutions.  Bawefveri  the  hlvarlate  gsna  and  blvaziate 
e.d.  that  have  been  suggested  have  been  produeed  by  the  saa»  teehnlgpe. 

This  teehnlque  also  produoes  the  bivariate  nomial  if  the  marginal  are 
noonal,  and  vlU  work  if  the  distribution  used  is  additive  (e.g.,  the  sum 
of  the  two  gasana  (or  e.d.)  randoax  variables  is  gasma  (or  e.d.)). 

Ibe  teehnlqoe  is  as  foUows.  Let 


“l  -  ♦  *3> 

(6t) -  - 

»J  ‘  ♦  *3^ 

Where  the  random  vaxlahles  and  3^  are  statistle^Jy  Independent 

and  have  the  appropriate  underlying  distribution,  l.e.  gaanw,  e.d.,  normal, 
ate.  Zf  the  dif^butlons.  of  Xj^,  ^  ewaw  (e.d.)  with  psraaaitsTS 

0^  and  (0^,82, a^)  respeetlvely,  then,  by  the  additive  pwpsrtlas 


32‘ 
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of  the 


i  gaaMi  (e.d.)  dlatrlhutlon,  the  distribution  of  vlll  be  sealed 
(e.d.)  vith  scale  factor  and  paraaster  (a^  -f  a^). 


StBllarlyf  the  distribution  of  Dj  vlll  be  sealed 


(e.d.)  vlth 


scale  factor  and  parameter  course,  and 

are  oozrelated  because  of  the  cobspu  random  variable  x^.  There  are 

five  parasHtars  involved,  ^3  these 

can  be  estlsmted  ftoai  and  ^  .  Such  a  procedure  Is 

1  J  1  J  I  i.  i 


described  In  [9]  ^br  the  o.d.  distribution  and  Is  the  basis  for  the 


bivariate 


used  In  tlO]. 


Bov  let  9^  be  the  underlying  distribution  of  ,  that  Is  let 


(68)  •j(t)  -  Rrobixj  <  tj 


vbare  I*  Is  the  cuaulatlve  gaaaa,  or  o.d.,  or  vhatever  ve  are  using. 


probsiblllty  distribution,  and  lot 


(69)  ^i(t) 


Bov  bgr 


(60)  sod  (^) 


53« 


(70) 


■  Prob(D^  <  <  Sg) 


•  ProlbtSj^^Cxj^  +  Xg)  <  Oj^j  SgCxg  +  Xj)  <  flgl 
Kov  the  ecndltlon  prohabllity  holdlne  x^  fixed  is 
(71)  Proh  [Sj^(x^  *3^  ^  *  *3^  -  ^2^*3^ 

-  Fl«b  tx^  <  -  X3,  *g  <  ^  -  *3  Ixjl 

0  0 

-  «l<87  -  »3> '^<1^  -  _ 

the  laet  step  hy  Bqiuitlon  (68)  and  the  IxMlepaiidenee  of  x^,  Xg,  and 


Hence  the  imoaodltlanal  prohSblUtp^  is 


-  *3>  3'^  -  *3’ 

. *3 

Itr  w  slailsr  WjfiMBt  we  tihtein 

(73)  >  e,;  Dj  >  eg) 

-  y'd  -^(Ji  -  t)]Cl  .  -  t)]  dir2(t) 

t  I  2 
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Differentiating  F  .  as  given  l)y  Equation  (72)  ve  obtain 

IjJ 

(7>.)  -  s;^  -  »>  ^  - »)  H'-'>  “ 

We  can  now  obtain  the  final  eiqpresslons  for  S(k^kj).  Substituting 
Equation  (74)  into  (Gv)  ve  obtain 

(75)  f  f  ■  t)((3(t)dtde^dS, 

®1  ®2  * 


Substituting  Equation  ('f^)  Into  (66)  >  ve  obtain 


(76)  E(k^kp  "Iff  - 1))  dw^(t)<u.(ej^)di.(eg) 

Jlk  ^  ^  2 


«1  ®2  " 
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